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Abstract 

We study the majority rule transformation applied to the Gibbs mea- 
sure for the 2-D Ising model at the critical point. The aim is to show that 
the renormalized hamiltonian is well defined in the sense that the renormal- 
ized measure is Gibbsian. We analyze the validity of Dobrushin-Shlosman 
Uniqueness (DSU) finite-size condition for the "constrained models" corre- 
sponding to different configurations of the "image" system. It is known that 
DSU implies, in our 2-D case, complete analyticity from which, as it has 
been recently shown by Haller and Kennedy, Gibbsianness follows. We in- 
troduce a Monte Carlo algorithm to compute an upper bound to Vasserstein 
distance (appearing in DSU) between finite volume Gibbs measures with dif- 
ferent boundary conditions. We get strong numerical evidence that indeed 
DSU condition is verified for a large enough volume V for all constrained 
models. 

Keywords: Majority-rule; Renormalization-group; non-Gibbsianness; Finite-size condi- 
tions; Complete analyticity; Ising model. 
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1. Introduction. 

In the recent few years many efforts have been devoted to the problem of a correct 
definition, on rigorous grounds, of various real-space renormalization-group maps. 
The main question is whether or not a measure 



arising from the application of a renormalization group transformation (RGT) T&, defined 
"on scale 6" , to the Gibbs measure \x is Gibbsian. In other words we ask ourselves whether 
or not v is the Gibbs measure corresponding to a finite-norm translationally invariant 
potential so that the "renormalized hamiltonian" is well defined. 

To be concrete let us suppose that n = /jp y h is the Gibbs measure describing 2-D 
Ising model at inverse temperature (3 and external magnetic field h ^ 0. Moreover we 
assume that our RGT can be expressed as: 



where T b (a' , a) is a normalized non-negative kernel. The system described in terms of the 
a variables by the original measure \i is called "object system". The cr"s are the "block 
variables" of the "image system" described by the renormalized measure v. 

We can think of the transformation T b as directly acting at infinite volume or we can 
consider a finite volume version and subsequently try to perform the thermodynamic limit. 
We refer to the basic reference [EFS] for a clear and complete description of the general 
set-up of renormalization maps from the point of view of rigorous statistical mechanics. 

The above mentioned pathological behaviour (non-Gibbsiannes of v) can be a con- 
sequence of the violation of a necessary condition for Gibbsiannes called quasi-locality 
(see [Ko], [EFS]). It is a continuity property of the finite volume conditional probabilities 
of v which, roughly speaking, says that they are almost independent of very far away 
conditioning spins. 

In [Ko] it is shown that a sort of converse statement holds true; namely: quasilo- 
cality + nonnulness (uniform positivity of conditional probabilities) of a stochastic field 
implies Gibbsianness but only in the sense of the existence of a finite norm but in general 
not translationally invariant potential associated to v. The construction of the poten- 
tial in Kozlov's proof is somehow artificial: it involves reordering of a semi-convergent 
sum. To get a translationally invariant finite norm potential one needs some additional 
stronger assumptions on how weakly the conditional probabilities of v depend on far away 
conditioning spins. We refer to [BMO] for a more detailed discussion on this point. 



v = T h \x 



(1.1) 




(1.2) 
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In some situations (see, for instance [C], [HK]) it is possible to use much stronger 
methods, based on cluster expansion, to compute renormalized potentials, showing finite- 
ness of their norm. 

In many interesting examples (see [El], [EFS], [EFK]) violation of quasi-locality and 
consequently non-Gibbsiannes of the renormalized measure v is a direct consequence of 
the appearance of a first order phase transition for the original (object) system described 
by fi conditioned to some particular configuration of the image system. More precisely 
given a block configuration a' let us consider the probability measure on the original spin 
variables given by: 



It defines the "constrained" model corresponding to a' (which here plays the role of an 
external parameter). 

For some particular a' it may happen that the corresponding measure [i a t (a) exhibits 
long range order. See also [GP], [I] where this mechanism was first pointed out. 

One can ask himself about "robustness" of the pathology of non-Gibbsianness. There 
are examples (see [M04]) in which, even though the measure v = Tbfip^h is not Gibbsian, 
one has that with the same (3, h, by choosing b' > b sufficiently large, the measure v' = 
Tb'/J>f3,h is Gibbsian. Alternatively one can think to iterate the map and, even though after 
the first step the resulting renormalized measure is not Gibbsian, it may happen, after a 
sufficiently large number of iterations that one gets back to the set of Gibbsian measures. 
This is often related to the fact that, given suitable values of the parameters f3, h (near the 
coexistence line h = 0, (3 > (3 C ), on a suitable scale b some constrained model can undergo 
a phase transition (somehow related to the phase transition of the object system); whereas 
given the same h, (3, for sufficiently large scale b any constrained model is in the one-phase 
region. 

Another notion of robustness of the pathology (see [LV] , [M05] ) refers to see whether 
or not it survives after application of a decimation transformation (see [EFS]); this can be 
relevant since decimation transformation does not change the thermodynamic functions or 
the long range correlations. 

Finally we want to say that weaker notions of Gibbsiannes of a renormalized measure 
v can also be considered. The usual notion of Gibbsiannes requires a control of quasi- 
locality of v uniform in a'. It may happen hat the particular a' responsible for the 
pathology is highly non-typical w.r.t. v. It appears plausible to ask for quasi-locality 
only for z/-almost all configurations a' (see [D2], [FP], [Ma]). We refer to [E2] for a nice 
up-to-date review of all the above problematic. 



Ha>(<T) 
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Now we want to notice that in many examples it happens that even though the 
object system is well inside the one-phase region, nonetheless for some particular block 
configuration a' the corresponding constrained model undergoes a first order phase transi- 
tion. Conversely there are many indications that if the constrained models are in the weak 
coupling regime then Gibbsianness of the renormalized measure follows. Recently Haller 
and Kennedy gave very interesting new rigorous results in this direction. They proved, 
under very general hypotheses, that if all constrained models are uniformly completely 
analytical (see [DS2], [DS3]) then the renormalized measure is Gibbsian with a finite norm 
potential which can be computed via a convergent cluster expansion. 

Let us give now the example of the Block-Averaging Transformation (BAT). Suppose 
to partition Z 2 into squared blocks Bi of side 2. In this case the new measure v is obtained 
by assigning to each block Bi an integer value rrii G {—4, —2, 0, +2, +4} and by computing 
the probability, w.r.t. the original Gibbs measure (ip^ of the event: J2 x eB, a x = m «- Then 
in this case we have: 



In this case a constrained model is a "multicanonical" Ising model namely an Ising model 
subject to the constraint of having, for every i, magnetization rrii in the block Bi. It has 
been shown in [EFK] that for BAT transformation the constrained model corresponding 
to rrii = V i, undergoes a first order phase transition at low enough temperature and 
that this implies violation of quasi-locality and then non-Gibbsianness of the renormalized 
measure v. Notice that for any constrained model with given {rrii} the value of the external 
magnetic field h is totally irrelevant. On the other hand, for h very large one can prove, 
by standard methods, absence of phase transition for the original model in the strongest 
possible sense: complete analyticity in the strong Dobrushin-Shlosman sense holds true 
in this case. This, as it has been shown in [EFS] gives an example of non-Gibbsiannes 
of a measure v arising from the application of a renormalization map to a measure \i 
corresponding to the very weak coupling region. We remark that, as it has been shown in 
[EFS], this non-Gibbsiannes is robust w.r.t. the choice of the scale b (or w.r.t iteration) 
whereas for large h it can be eliminated by applying one decimation transformation (see 



Other interesting examples have been found (see [EFK], [El], [E2]). 

We want to stress that, in general, it is not sufficient to control that one single 
constrained model is in the one-phase region to imply Gibbsianness of the renormalized 
measure. In [CG] and in [BMO] for the BAT transformation it was suggested that the 




[M05]). 
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fact that the constrained model with nii = was in the high temperature phase could 
be sufficient to imply the existence of a finite norm renormalized potential. Recently A. 
van Enter showed with an example that this believe is not sufficiently justified and some 
extra arguments related to the specific nature of BAT transformation are needed to imply 
Gibbsiannes from absence of phase transitions for {nii = 0} constrained model (see [El], 
[E2]). 

Thus, in general, the moral is that what is relevant for Gibbsiannes of v are the 
(intermediate) constrained models; we repeat that it can be sufficient that even only one 
constrained model undergoes a phase transition with long-range order (despite of the 
possible very weak coupling regime of the object system) to imply non-Gibbsiannes of v\ 
whereas, in general, absence of phase transition in a very strong sense is needed for all 
constrained models to imply Gibbsiannes of v in the strong, cluster expansion, sense as it 
has been shown in [HK]. 

In the present paper we will analyze a particular RGT: the majority rule transforma- 
tion applied to the 2D critical Ising model. It will be precisely defined in next Section 2. 
This transformation, in the same situation of criticality, has been studied by Tom Kennedy 
in [K2]. The author establishes some rigorous results reducing absence of phase transi- 
tion for some particularly relevant constrained models to the verification of some suitable 
"finite size conditions" introduced in [Kl]. It is an "almost computer assisted proof of 
absence of phase transition for these constrained models. 

In the present paper we inquire for the validity, in principle for every possible con- 
strained model, of a finite size condition: the Dobrushin-Shlosman uniqueness condition 
(DSU, see Section 2 below) which implies (as it will be explained in Section 2) complete 
analyticity and then, using the results of [HK], Gibbsianness of v. 

Strictly speaking the proof of [HK] does not directly apply to our case since it requires 
the condition that the kernel T(a',a) is strictly positive for every a', a. Probably this is 
only a technical restriction that can be removed ([K3]). In any case in [HK] the authors 
claim that for the majority rule they are able to obtain an equivalent system with T(a' , a) > 
by first summing out some spins in the original system. 

Our results and their strength are, in a sense, complementary to the ones of [K2]. 
Our study will be numerical but, similarly to [K2] not only in the sense of "traditional" 
Monte Carlo simulations. Rather, for each constrained model, we will try to measure by 
a computer a quantity appearing in DSU such that if we could rigorously prove that it is 
strictly less than one then we could deduce from some theorems strong properties typical 
of the one-phase region, for arbitrarily large and even infinite systems. 

In [K2] some constrained models were analyzed in terms of a finite size condition 
easier to be satisfied than DSU, for which the author could also have provided a compu- 
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tation based on interval arithmetics suited for a computer assisted proof. This finite size 
condition of reference [Kl] is not sufficient to imply complete analyticity. In the present 
paper, as we said before, we try to verify for every constrained model DSU condition; 
but, as it will appear clear in Section 5 we can, with nowadays machines, only perform a 
Monte Carlo calculation. For many reasons we cannot, at the moment, hope to improve 
our calculations to get a complete control and possibly a computer assisted proof. We 
will explain in Section 6, devoted to the conclusions, in what sense our results can be 
considered satisfactory. 

The DSU condition involves the calculation of the so-called Vasserstein distance 
between two Gibbs measures in a finite volume with boundary conditions differing only 
in one conditioning site. In a recent paper ([BMO]) the authors, in the context of BAT 
transformation, for one particularly relevant constrained model ({to;} = V i) tried to 
verify the same DSU condition but they were only able to provide a (numerical) lower 
bound for the concerned Vasserstein distance. The reason was that a lower bound (see 
(2.16)) in terms of total variation distance can be found involving thermal averages; thus 
this lower bound is well adapted to be studied by Monte Carlo methods but, on the other 
hand, it is only able to give some indications on the validity of the true condition and this 
since it appears reasonable to expect that it is a good lower estimate; strictly speaking it 
is only useful to disprove the condition. 

In the present paper we present a Monte Carlo algorithm, inspired by the "surgery 
method" introduced by Dobrushin and Shlosman (see [DS1], [DS2]) that we call dynamical 
surgery. It provides an upper bound to the Vasserstein distance and so it goes into the 
correct direction to prove DSU condition. 

Our numerical results strongly suggest that indeed DSU condition is satisfied in our 
present situation for all constrained models and this, as we said before implies Gibbsiannes 
of renormalized measure. 

The paper is organized as follows: in Section 2 we define in detail our majority rule 
transformation and the constrained models. In Section 3 by using "conventional" Monte 
Carlo methods we provide rough estimates of the critical temperatures of some particularly 
relevant constrained models. In Section 4 we introduce our algorithm. In Section 5 we 
give our main numerical results. In Section 6 we give the conclusions. In Appendix A we 
present the computation of the best joint representation of two measures which is used in 
our algorithm. 
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2. The majority rule transformation and the constrained models. 

We will consider the usual (ferromagnetic nearest neighbors interaction) 2D Ising model, 
with zero external field, on a finite square with even side: A := {1,...,2L} 2 C Z 2 with 
LeW. We denote by a G Oa := { — 1, +1} A a configuration of the system in A and by 
r G { — 1, +1} 9A+ a boundary condition, namely a configuration in <9A+ defined by 

<9A + := {i G Z 2 \ A : 3j G A : % and j are nearest neighbors} . (2.1) 

By ai G {— 1, +1} and tj G {—1, +1} we denote the spin variables on the site i G A and 
j G dA + . It is also convenient to think of r as an extended configuration in { — 1, +1} Z \ A . 

The energy associated to a G Oa with r boundary condition outside A and zero 
external magnetic field is given by 

h t a {<j) •■=- - a ^ w G ° a ' ( 2 - 2 ) 

<i,j> <i,3> 
iJeA ii£ A,je8A+ 

where the first sum runs over all pairs of nearest neighbors sites in A, while the second 
sum runs over all pairs < i,j > of nearest neighbors sites such that i G A and j G <9A+. 
The Gibbs measure describing the equilibrium properties of the system at the inverse 
temperature (3 is denoted by fx^ A (a) Va G Oa and is given by 

e -f3Hl(a) 

VpA*) : = ^ p -(3HUr,) VffG °A- ( 2 - 3 ) 

In the following four steps we give the precise definition of the "Majority Rule" (on 
scale 2) transformation. 

1. Vx, y G Z let denote the 2x2 block whose center has coordinates: 

(2x — |, 2y — -|); the collection of all blocks -B( x , y ) Vx, y G Z gives rise to a partition of the 
lattice Z 2 . If we restrict ourselves to pairs (x, y) G {1, L} 2 we get a partition of our box 
A. Given a G Oa we denote by &\ x y y •••i (T f X y) the four spins corresponding to the four 

sites of the block -B( x , y ) and we define m^ x ^ := Y^t=i a \x y y we su PP ose the four spins 
a {x y) ^ = ordered in lexicographic way. 

2. We define the new lattice A' by collecting the centers of all -B( x , y ) blocks and by rescaling 
the lattice spacing by a factor two; the site of A', which is the center of the block B^ x y ^ 
will be simply denoted by the pair (x, y). 
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3. On each site (x,y) G A' we define the renormalized spin o~^ x ^ G { — 1, +1} and we 

consider the space n' A , := {-1,+1} A \ We define the kernel K : (a, a') G A x Q' A , -> 
K(a,a') G {0,1} as follows 

if 3 (x, y) G A' : m (Xiy) 7^ and m( x , y ) • ^J^) < 
K(a, a') := \ if 3 (*, y) G A' : m (x , y) = and a^ y) ■ a'^ y) < . (2.4) 

1 otherwise 

4. The Majority Rule Transformation is the transformation which maps the "object" 
model (A,£l\, ftp A (a)) onto the "image" model (A', Q' A ,, /j,'p A ,(cr')), where 

^ A ' (a) '= E^E^ A ^,0^ A (.) WG °-- (2 - 5) 

Notice that we could have used the notation Q A i in place of Q' A , since here, for the majority 
rule transformation, contrary to other transformations like BAT, the single renormalized 
spin variable still takes values in {—1, +1}. 

This transformation is well known in physics literature and it has been widely used 
to investigate the properties of many spin models (see, e. g., [NL] and references therein). 

A very important role in our discussion will be played by the "constrained models" : 
given a' G Q' A , we call constrained model corresponding to a', and we denote it by Zp a i, 
the model (A, Q A , yTp a a '( a ))' where 

^ Aa) = E* 0t (2 ' 6) 

is a probability measure on Vt A . Due to the fact that K(a,a') G {0,1} we have that, 
Vcr G Oa and Vcr' G A , the constrained model Tp a , can be seen as the model 

A.V^W'W = ^ e -Pm( v ) > ( 2 - 7 ) 

where we have introduced the constrained configuration space 

n A ,a' := W G A : V(x,y) G A' : 
(m (a ., y) ^ and m (a!)V ) • (7^ ^ > 0) or (m (a . )V) = and a\ x y) ■ a[ x y) > 0)} . (2.8) 
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In other words we can say that the constrained model Xg a , is a model defined on the origi- 
nal lattice A, with the same hamiltonian HJJ^a) as the object model, but with configuration 
space Oa,<t'. 

The microscopic states of the constrained models can be characterized by means of a 
suitable block-variable. In the original Ising model there are 2 = 16 allowed configurations 
in each block B( x y y We partition these block-configurations into two disjoint classes C + 
and C- as follows: 

• block-configurations belonging to class C + 



- + 




+ - 




+ - 




+ - 




+ + 




+ + 




+ + 




+ 


+ ' 


+ + 




- + 




+ - 




+ + 








- + 




+ - 




+ 


+ _ 



(2.9) 



• block-configurations belonging to class C_ 



















- + 




- + 




- + 




+ - 






- + 




+ - 




+ + 








- + 




+ - 







(2.10) 



For any constrained model 2^ a , we allow, in the block B^ x y y only the configurations in 
the class C s ^„ na i ; in order to classify these block-configurations we introduce the block- 
variable S^ x ^ G {1,2,..., 8}: to the values S^ x ^ = 1,...,8 correspond, respectively, the 
eight block-configurations in (2.9) if a'^ x ^ = +1, the eight block-configurations in (2.10) 

if = - L 
Warning. 

Here, and in what follows, we use the (non-conventional) expression block-variable referring 
to the variable S^ x ^ taking values in the set {1, 2, 8}; the spin variables a'^ x y y defined 
on the renormalized lattice A' and taking values in {— 1,+1}, will be sometimes called 
renormalized variables. 

Obviously, given any constrained model Xg a , , one and only one configuration S G 

Oa' := {1,2, ...,8} A can be associated to any a G fiA )(T ', and viceversa. Hence, each state 
of the model Xg a , can be represented by a collection S G 0/y of block- variables S( XjV ), but 
we recall that the "meaning" (in terms of the original spin variables) of each block-variable 
S(x, y ) depends on the sign of <j'^ x y y 

Due to the bijection between Oa,<t' an d the hamiltonian HJ^(a) of the model 
Xg a , can be thought, for any a G Oa,<j', as a function B.\, a , T ,(S) of the block- variables 
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configuration S G Oa'! here by r' we mean a configuration of the renormalized variables 
in dA' + , that is r' G fi' 9A ,+ := {-1,+1} 9A ' + , where 

dA ,+ :={(x,y)eZ 2 :0<x,y<L+l, x ^ y, (x,y) ? A'} ; (2.11) 

by T G Qqa>+ '■= {1) ...,8} 9A + we mean the boundary condition expressed in terms of the 
block-variables in the set dA' + . Notice that given the original r and r' in general T is not 
uniquely determined. 

We remark that a block- variable in Br Xjy \, for (x,y) G A' U <9A /+ , is completely 
"meaningless" if the correspondent renormalized variable has not been specified. Finally, 
one can say that the constrained model (2.7) is equivalent to the model defined on the 
lattice A', with configuration space Oa' and equilibrium measure given by 

-PH*, , ,(S) 

»b.a:*>.AS)--= m V ^ G °A' • (2.12) 

Now, we want to study the above defined constrained models by means of the finite 
size Dobrushin-Shlosman uniqueness condition (DSU); to introduce DSU we need some 
definitions. 

Let us consider two measures pi and p 2 on a finite set Y; let p(-, •) be a metrics on 
Y and denote by fC(pi, p 2 ) the set of joint representations of p\ and p 2 , namely the set of 
all measures p on the Cartesian product Y xY such that 

p(y,y') = Pi(B) and Ky,y') = V2(B) vbcf . 

y eB,y'eY y eY,y'eB 

We set: 

Var{p u p 2 ) := 2 l^ 1 ^ ~ ^ 2 ^)l ( 2 - 13 ) 
yeY 

and, given a metrics ponY 

v p (p u p 2 )-.= inf ^ p{y,y')-»{y,y') ; (2.14) 

M6/C(Mi,M2)^, 6y 

Var(//i, /X2) and V p (pi, p 2 ) are respectively called iota/ variation distance and Vasserstein 
distance with respect to p between the two measures p\ and p 2 . 

Let us consider a spin system on Z d with single spin space S and range-one interac- 
tion (generalization to arbitrary finite range is trivial); we denote by rji the spin variable 
associated to the site i G Z d . 
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For any finite set V C Z d we denote by dV + the set of points outside V whose spins 
interact with the spins inside V and by r\ v G <S y a spin configuration on V. Given the 
boundary condition £ G iS v , we denote by fj,y the Gibbs measure in V with boundary 
conditions £ outside V. Given a metrics p on 5 we associate to it a metrics p v on <S^ 
defined as follows 

PviVv, Vv) ■= P^n V'i) VVv, Vv ■ 
iev 

We say that condition DSU P (V, S) (the Dobrushin-Shlosman Uniqueness condition 
in V with respect to the metrics p and uniqueness parameter S) is satisfied if and only if 
3 a finite set V C Z d and 35 > such that: Vj G <9y + 3a j > such that for any couple 
of boundary conditions £, £' G iS^ with ^ = £• Vz 7^ j one has 

and 

Let us consider the metrics on the single site variable given by: 

,1 iff 77 7^ 77' 

#77,7/) :=<j /7 ^/ ; (2.15) 

U otherwise 

in this case condition DSUp(V, 5) will be simply denoted by DSU(V, S). 
We observe that, in this case 

2vo4,/4) > Y. V M,4) = Y, Va M>4) ( 2 - 16 ) 

iev iev 

as it easily follows from (2.14) and Proposition Al in the Appendix. 

In order to describe Dobrushin-Shlosman' s results based on DSU p (V,S) condition 
we need some more definitions. We are going to introduce two kinds of mixing conditions, 
of a priori different strength, for Gibbs measures in a finite volume A. The first ones are 
weak mixing (WM) conditions saying that the influence of a local change in conditioning 
spins decays exponentially fast with the distance from the boundary <9A; the second kind, 
the one of strong mixing (SM) conditions, corresponds to the case where the influence 
of a change of a conditioning spin x decays exponentially fast with the distance from x. 
We refer to [MOl], [M02] for a critical discussion of these different notions. It has been 
shown with some examples (see [Sh]) that it may happen that WM is satisfied whereas the 
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corresponding SM is not. Since we are speaking of finite volume mixing condition we have 
to make explicit in the notation the constants referring to the concerned exponential decay. 
Of course a particularly interesting case is when we have these mixing conditions in a class 
of arbitrarily large volumes A with uniform constants. We refer again to [MOl], [M02] for 
a discussion of these point. It turns out that a crucial aspect is the class of volumes that we 
are considering. In particular in the Dobrushin-Shlosman theory of complete analyticity 
(see [DS2], [DS3]) arbitrary shapes were considered whereas in the approach developed in 
[O], [OP], [M02] only sufficiently regular domains were involved. 

We say that a Gibbs measure on fi A satisfies a strong mixing condition with 
constants C, 7 if for every subset Ac A: 

sup Var(^ AA ^a) < Ce~^ dist ^ (2.17) 

where Tx = r x for x 7^ y and fx A A is the /x A -probability distribution of the spins in A. 
We denote this condition by SM(A, C, 7). 

We say that a Gibbs measure /i A satisfies a weak mixing condition with constants 
C, 7 if for every subset A C A 



sup Var(/j, T AA , nl A ) < C ^ exp(-7|x - y\) ; (2.18) 

x£A, yedA+ 



r,T'en A 



we denote this condition by WM(A, C, 7). 
Theorem DS (see [DS1]) 

Let DSU ( V, S) be satisfied for some V and 5 < 1; then 3 C > 0, 7 > such that condition 
WM(A, C, 7) holds for every A. 

Theorem MOS (see [MOS]) 

Let the dimension of the lattice be d = 2. If there exist positive constants C and 7 such 
that the Gibbs measure /i A satisfies the weak mixing condition WM(A, C, 7) for any finite 
A C Z, then there exist positive constants C and 7' such that the Gibbs measure /x A 
satisfies the strong mixing condition S'M(A,C",7 / ) for any sufficiently regular domain A 
and in particular for any square A^ with arbitrary side L. 

Here "sufficiently regular" means "multiple of a sufficiently large square" (see [MOl]). 
Remark 

Actually the conclusion of the above theorem remains true even if we assume the weak 
mixing not for all finite subsets of Z but only for all subsets of a square A^ o provided that 
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L is large enough (depending of course on the constants C and 7 and on the range of the 
interaction) . 

Then our strategy to show Gibbsianness of the renormalized measure arising from 
the application of the Majority Rule transformation is based on the following chain of 
implications. 

We try to verify DSU(V, S) for some given V and 6 < 1 for any constrained model. 
Then if we could apply Theorem DS we would get WM( A, C, 7) VA with the same constants 
C, 7 for all A and for all constrained model. Subsequently, if we could apply Theorem MOS 
(by exploiting two-dimensionality), we would get SM(A, C, 7') for all sufficiently regular 
domains A with the same constants C',7' for all these A and for all constrained model. 
This would directly imply the validity of the conclusions of Theorem 1.1 in [HK]; indeed 
it immediately follows from the proof of Theorem 1.1 in [HK] that the authors could have 
obtained exactly the same result by only assuming their strong mixing hypothesis, (that 
they express in an equivalent form, valid for Ising-like systems, as exponential decay of 
two-points spin-spin truncated correlations) uniformly in the constrained model and in the 
boundary conditions only for all sufficiently regular (in the above specified sense) domains 
instead of for all domains. 

Strictly speaking Theorems DS and MOS apply to translationally invariant situations 
and our constrained models are not, in general, translationally invariant. However it 
is easy to convince oneself that both Theorems extend in a straightforward way to the 
non-translationally invariant case provided we assume spatial uniformity of the bounds 
appearing in the hypotheses. In other words we have that it is sufficient to assume the 
validity, for a given constrained model, of DSU(V, 5) for some V and some S < 1 uniformly 
in the location of V to imply SM (for arbitrary sufficiently regular volumes with uniform 
constants) for the same model. But if we are able to show that DSU(V, 5) is verified for 
a given V and S < 1 for all constrained models (namely for all a' G fiyi) this implies, at 
the same time, (using extended versions of Theorem DS and MOS) the validity of SM (for 
arbitrary sufficiently regular volumes with uniform constant) for all constrained models 
and then, via Theorem 1.1 in [HK], Gibbsianness of the renormalized measure. 

Now we start applying the above strategy. We first introduce some specific defi- 
nitions. Let us consider a squared volume V C A' with side / and the corresponding 
subset V, with side 21, of the original lattice A. Let us consider the metrics on the single 
block- variable space {1, ...,8} defined by 





otherwise 




(2.19) 



1/ febbraio/2008 [13] 



2:7 



and the metrics on fly' given by 

Pv ,(S,S'):= K s i*,v)> s {x,„)) VS,S'en v >. (2.20) 

(x, y )ev> 

Now, given a' G fly, and r' G Q' dv , + , we denote by [Ti, T2]( x ,y) a pair of boundary 
conditions 7i, T 2 G Q 9V ,+ such that 7i (x , V) = T 2(x , yl) V(x',y') G <9F /+ \ {(x,y)} and we 
define 

c-[7i,T 2 ] (Xiy) , Tl Ta |dF' + | _ Tl Ta 4 

(2.21) 

where , (a*/jV' <t' t'' ^/?V' <t' t') * s ^ ne Vasserstein distance between two equilibrium mea- 
sures for the constrained model corresponding to a', t' which have been obtained by mod- 
ifying the boundary conditions just in one site in dV' + . We set: 

Spy, a' ■= sup sup EgyT^J^ ; (2.22) 

(x,y)edV'+ [Ti ,T 2 ]( Xt y) ' ' ' 

heuristically we can say that Spy j(7 i measures how much the equilibrium of the constrained 
model is modified if one changes the boundary condition in one site, uniformly in the site 
and in the boundary conditions; it can be called uniqueness parameter (in the sense of DSU 
condition) in V for the constrained model characterized by a' at inverse temperature (3 

It is a trivial consequence of the definitions that if Spy jCr > < 1 then the Dobrushin- 
Shlosman uniqueness condition DSU(V, 6) is satisfied for some 5 < 1 for the constrained 
model corresponding to <t',t'. The main aim of this paper is to build up a Monte Carlo 
algorithm in order to estimate Spy , a ' and to show that at (5 = (3 C : = | log(l + y/2) (the 
critical inverse temperature of the standard 2D Ising model) there exists a volume V such 
that for all possible constrained models Spy ^> < 1, namely: 

Spy := sup £p,v', a > < 1 • (2.23) 
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3. Rough estimate of critical temperatures. 

In the following sections we will see that two constrained models (the chessboard and the 
striped model, see definition below) are particularly "dangerous" , that is they satisfy DSU 
condition for volumes larger than those needed by the other constrained models. 

The chessboard and the striped models are the constrained models corresponding, 
respectively, to the two configurations C,5 6 fi' defined as follows 



these two configurations are respectively depicted in Fig. lb and Fig.lc. 

We study the equilibrium properties of the chessboard and the striped model by 
means of a standard Monte Carlo procedure based on a suitable Heat Bath dynamics. 
This dynamics is generally used to compute mean values (with respect to the equilibrium 
Gibbs measure) of some observables as time averages; this mean value will be denoted by 
< • >. 

This "conventional Monte Carlo" analysis is preliminar to the main numerical results 
of this paper. Here we are not interested in a large scale simulation nor in a precise estimate 
of critical points and exponents; we just want to have a strong evidence of the fact that the 
critical inverse temperature of the two dangerous models is, in both cases, much greater 
than [3 C . 

By making use of the notation introduced in Section 2 (relatively to the constrained 
models Xg a ,) we describe, now, the discrete time Heat Bath dynamics used in our Monte 
Carlo study. It is given by the Markov chain defined below: 

1. We consider the constrained model corresponding to a' G A , and we assume periodic 
boundary conditions; we denote it by I/3, a ' and its hamiltonian by Ha> j(7 '(S) with S G Qa'- 

2. At each step we perform a complete updating of all S( x>y ) block-variables following 
the lexicographic order; in the remaining of this section this will be called a Monte Carlo 
sweep. 

3. For any (x, y) G A' the new value S', x s of the block- variable S( x>y ) is chosen at 
random according to the Gibbs measure in -B( x , y ) with boundary condition S^ x y ^ := 
S\(A'udA' )\b, n- Hence, if we denote by H a >(Si „\\S? s) the contribution of the block 
B(x,y) to the energy of the system, the transition probability is given by 




(3.1) 



(3.2) 



E 8 , 
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We remark that in our notation there is'nt any explicit dependence of the boundary con- 
ditions (T-t'), because they are supposed to be periodic. 

In order to estabilish for which value of the inverse temperature (3 the two models 
Xg ; c and Tf},s are critical, we have computed the specific heat defined in terms of the 
equilibrium energy fluctuations: 

C A ,<r> ■■= j^(< Hl ja , >-< H A ,, a . > 2 ) W = C,S . (3.3) 

Both models have been studied in the case L = 64, which means that we have 
considered the two models defined on square lattice containing 128 2 original sites; we have 
considered smaller values of L, as well, in order to check the finite size behaviour. In the 
case of model Zp,c we have performed 10 5 full sweeps of our Monte Carlo algorithm for 
each value of /3, while in the case of model 1.5 x 10 5 sweeps have been performed. 

In Fig. 2 and in Fig. 3 we have plotted the specific heat as a function of the inverse 
temperature (3 in the case of models Ip,c an d 1p,s respectively; these results have been 
obtained by analyzing the Monte Carlo data by means of the "jackknife" procedure. From 
the pictures it is clear that both inverse critical temperatures, (3^ and (3f , are significantly 
greater than the Ising inverse critical temperature (3 C . Our rough estimates are the following 
ones: 

$ = 1.60 ±0.05 /3f > 2.2 . (3.4) 
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4. The algorithm. 

Let us consider the volume V introduced in Section 2, a' G fly/, r' G fl' dv , + , a pair of 
boundary conditions [Ti,T 2 ]^,y), the two equilibrium measures fJ^ A , a , T , and fJ^ 2 A , CT / T ,- 

Let us denote by K^J V T 2 ^ X ^ ) the set of joint representations of AtJ 1 ^ a , T , and fJ^ 2 A , a , r ,\ 

we want to give a numerical estimate of the quantity E^]ff 2 a ^^ defined in (2.21). 

A similar problem has been studied in [BMO]: in that paper, by means of a standard 
Heat Bath dynamics, a lower bound to £g ^, ^./^^ has been calculated; now, we build up 
a Monte Carlo algorithm in order to obtain a numerical estimate of an upper bound to 
the same quantity. 

First of all we heuristically describe the idea on which the algorithm is based. In 
order to exactly calculate the Vasserstein distance one should know a joint probability 
measure Q*(-, •) G /Cg v , 2 J X T V on the space fly' x sucn that 

^ v >£vV^A*?V W ) = Q*(^ (1) ^ (2) )P^(^ (1) ^ (2) ) , (4-1) 

s( 1 ),s( 2 )en w 

that is a joint probability measure "optimizing" the sum in (2.14); this is a very hard 
task. On the other hand it is possible to calculate the Vasserstein distance between two 
equilibrium measures, relative to two different local boundary conditions, of a single block- 
variable; indeed, in this case, from the definitions given in Section 2, it turns out that 
the distance between two configurations is p(a, a') = 1 — 5 a , a i Va, a' G {1,...,8} (see 
(2.19)) and so the Vasserstein distance coincides with the total variation distance (see [Dl], 
page 472). In this case it is also possible to calculate the "optimizing" joint distribution 
measure. Hence the idea is to build up a dynamics describing the evolution of two coupled 
systems, with different boundary conditions, such that at each step a single block- variable 
is updated according to the "local optimizing" joint distribution law. 
Now we give the detailed definition of the dynamics. 

1. Given a' G 0,' v ,, r' G Q' dv ,+ and a pair of boundary conditions [7i,72](s,y), we want 
to describe the evolution of the pair of copies of the constrained model corresponding to 
a' G fly, and r' G {—l,+l} dv with boundary conditions respectively given by T\ and 

r 2 . 

2. We consider the Markov Chain X t := (S^(t), S™(t)) Vt G N where S^(t) and S^(t) 
are the configurations of the two copies of the system at the time t. The two processes 
S^\t) and S^ 2 \t) are called marginal chains, while X t is called joint chain. 

3. At each instant t one and only one site (x, y) G V is taken into account; all sites 
(x,y) G V are visited in lexicographic order, hence in an interval of time At = \V'\ all 
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sites {x, y) G V are considered: this is called a sweep (or, sometimes, full sweep) of the 
dynamics. 

4. Let us consider the site (x, y) which is updated at the time t; we introduce the following 
notation: Q Ti [S(J y )\Syi\{( x y)}^ — 1)] with % = 1,2 is the equilibrium Gibbs probability 
measure for the block- variable in (x, y) conditioned to the values of the block-variables in 
V \ {(x, y)} at time t — 1; we denote by /C( x y )(t) the set of joint representations of 

9 T '[ S S»)l S v.\ (( „ rt) ('-l)] and ^lS^\ lM} (t-l)) (4.2) 
and, finally, by 

q * Tl > T *[S ( l] y) ,S^ y) \Sy)\ {{Xty)} (t - 1)), Sy\ {(x y)} (t - 1)] (4.3) 
the joint representation in K^ x ^{t) such that 

^(« T, iC)i s< A««)>(«- 1 )]'9 T, [ s S.)i s ^«(«)>(«- 1 )]) = 

8 

P^(x,y)^(x,y)) 1 l D (x,y)^(x,y)\ D V'\{(x,y)}y l L )) ' °V'\{(x,y)} V ^J> 

s (l) ^(2) =1 
O.I/) 

(4.4) 

the joint representation (4.3) can be exactly calculated as it is shown in Appendix A. 
All probability measures introduced above depend on f3,V',a' and r'; we drop explicit 
dependence to simplify the notation. Notice that, of course, the q Ti [■{■] are stationary. The 
dependence on t — 1 in the conditioning spins is made explicit only to clarify the updating 
rule from the configuration at time t — 1 to the one at time t. 

5. Given t > 1 and the corresponding site (x,y) G V, the configurations S^(t) and 
S( 2 \t) are obtained by choosing the pair (S^ y ^(t), S^ y ^(t)) at random according to the 

joint probability measure (4.3) and S^, y ,^(t) = S^, yl ^(t — 1) Vi = 1,2 and V(x',y') G 
V'\{(x,y)}. 

In analogy with the terminology used in [DS1], we call this dynamics Dynamical 
Surgery; now some easy properties of this dynamics are discussed. 

1. Given t G N, it is immediate to see that if 

^S' ,!/')(* ~ X ) = ^S'.w') (* ~~ X ) V ( x '' y') nearest neighbors of (x, y) 
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then 

in other words this means that in this case the joint representation (4.3) lies on the diagonal. 

2. The evolution of the two marginal chains S^(t) and S^ 2 \t) is described by a stan- 
dard Heat Bath dynamics; this is a consequence of the fact that the probability (4.3) 
is a joint representation of the single site probability distributions (4.2), which are the 
Gibbs measures of a single block-variable with boundary conditions given, respectively, by 
{^i^v^iix^)}^ _ !)) and (^2, Sv'MixM)}^ ~ Hence tne equilibrium distributions 

for the two mare; inal chains S^\t) and S^ 2 \t) are the Gibbs measures with boundary 
conditions given, respectively, by T\ and T 2 . 

3. Given t G N, given Q(-, •) G K,f^T^;f\ if the stochastic variable X t = (S^ (t) , S<® (t)) 
is such that 

P(S^(t) = S, S<- 2 \t) = S') = Q(S, S') V5, 5' G JV , 



that is it is distributed according to the joint probability measure Q(-, •), then it is easy 
to prove that X t+ i is distributed according to an element of Kp A ', 2 J X T J\ as well. 



From the fact that the joint chain X t is an ergodic, aperiodic chain with a finite 
state space and from properties 2 and 3, it immediately follows that there exists a unique 
equilibrium distribution for the Markov chain X t ; it is a joint probability measure Q(-, •) G 
^■p A' ct' t'' - This means that, provided the two marginal chains have thermalized, that 
is they are almost described by the corresponding Gibbs measures, the joint chain, after 
a long enough time, is almost characterized by the distribution Q(-,-), which is a joint 
representation of the Gibbs measures Hp^, a , T >(S) and Hp 2 A , a , T >(S). Hence, if one needs 
to calculate the "phase average" of an observable depending on S^ 1 -* and with the joint 
probability Q, one can perform a sufficiently long Monte Carlo run with the dynamical 
surgery and calculate the "time average" of the same observable. 

If it were Q = Q* then we could perform a Monte Carlo calculation of the uniqueness 
parameter Sp y, , but in general Q ^ Q* and moreover we cannot say anything on 
how much Q is close to Q* . Anyway we can use the joint representation Q to calculate an 
"upper" estimator of the uniqueness parameter ; we set: 

uf: v T^ T V := \ - yl>< (S W (S- \V'\),S( 2 \s. \V'\)) (4.6) 

s=l 

where / is the number of sweep performed during the Monte Carlo run. It is clear that 
Up T yJ' 2 ( j^ T f ) is a numerical upper bound to sjy^jt*^ . 
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The same Monte Carlo run can be used to estimate a "lower" bound to the uniqueness 
parameter like in [BMO]. Indeed, we define the quantities 



1 i 

: = jE^c^^ivi) V * = 1.2V(x,y) EV'Vae {1,...,8} (4.7) 

s = l 



and we set 



(sB,»)ev ae{i,...,8} 
see [BMO] for more details. 
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5. Upper and lower bound: numerical results. 

In Section 4 we have developed a Monte Carlo algorithm in order to evaluate an upper 
and a lower bound on Ep y, 2 J X T J ] and we have denoted these two quantities respectively 

D y U /3,V',<t',t> and L '(5,V>,a>,T> • 

This algorithm cannot be used to give a direct evaluation of an upper bound to Epy\ 
indeed, in order to calculate an upper bound to Epy we should consider all possible con- 
strained models and all possible pairs of boundary conditions [71, T^\{ x ,y) ( see the notation 
introduced in Section 2), that is we should perform 3 • / ■ 2 l + 8l+1 runs of our Monte Carlo; 
one immediately realizes that this is an almost impossible task since it would give rise to 
an enormous computation even in the case of small volumes (I = 2, 3). 

Indeed, two boundary conditions can differ in one site (x, y) G dV' + in six ways: if 
one considers the block -B( x , y ) with (x, y) G dV' + and gives the block-variable S( XiV ), only 
two of the original spins a^ x y y those adjacent from the exterior to the lattice V, influence 
the equilibrium properties of the system; hence the possible ways in which two boundaries 
may differ are the following ones 

++ +- 
++ -+ 

+- — 
-+ — 

where we have depicted the two couples of original spins, belonging to the block -B( x , y ) and 
adjacent from the exterior to V, that one can have when one considers two boundary con- 
ditions differing just in (x, y). Hence, one can easily convince himself that the total number 
of pairs of boundary conditions [Ti, T^u.^ is given by 6-2 2 ( 4I-1 >; finally, noticed that given 
I there exists 2 l possible constrained models and that the site (x, y) can be chosen in 41 
different ways, it is immediate to see that considering all possible constrained models and 
all possible pairs of boundary conditions amounts to examine 3 • I ■ 2 l +8l+1 different situa- 
tions. This number could be reduced by taking into account various symmetries; however 
one would still have, for interesting values of /, an excessively big computation to perform. 

Hence we are forced to consider just some of the possible constrained models and some 
of the possible boundary conditions, that is we are forced to perform a sort of "statistics" . 

YT T 1 

Let us denote by Upy the highest value obtained for y, ^J^V and by Cpy the 

corresponding estimate of cJ^^jT^? . 

First of all we made a preliminary evaluation of our upper and lower estimators for 
different values of I, by choosing completely at random a certain number of constrained 
models and some of all possible pairs of boundary conditions. 
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We describe, now, how the statistics has been performed in all cases; our numerical 
results are summarized in Table 1. 

• 1=1: In this case the volume V contains a single site, that is V = {(1,1)}. The 
dynamics defined in Section 4 is based on the local updating of a single block-variable; 
this updating is worked out according to the probability distribution (4.3). Hence, in the 
case / = 1 the estimator (2.23) has been exactly calculated as follows: 

Spy. = 4 ( SUp SUp SUp J2 1* ^ l S (ll) > 5 (1 2) 1)] ) ; 

V'en' v , (x,y)eav+ [t u t^ v) (1) {2) / 

(i,i)' (i,i) 1 ""' ' 

q(l) -Cj( 2 ) 

s (i,ir (i,i) 

On the other hand since in this case the Vasserstein distance coincides with the total 
variation distance (see Appendix A) we could compute the same quantity by means of the 
expression (2.13). We got: 

£p,v = 2.119 

this result means that in this case there exists at least one constrained model (we remark 
that there are only two constrained models, respectively corresponding to a'^ ^ = +1 and 
a'^ ^ = —1) and one pair of boundary conditions [Ti,^]^^) such that the Dobrushin- 
Shlosman condition is not fulfilled. Hence, the volume V with I = 1 is not "large enough" 
for our purposes. 

In both cases cr^ ^ = +1 and cr^ ^ = —1 and for many pairs of boundary conditions 
\T\ , T 2 ] ( x , we have evaluated the estimator (2.21) by means of the Monte Carlo algo- 
rithm, as well; in this way we checked our computational procedure to get the best joint 
representation (4.3). The results that we obtained differ from the exact values by 1-2%; 
this shows that, at least in the case I = 1, our Monte Carlo procedure is very efficient. 

• 1 = 2: we have considered all possible constrained models and for each model we have 
considered 200 different pairs of boundary conditions. By performing 1.3 • 10 6 full sweeps 
of our Monte Carlo, we have obtained the results in Table 1; these results show that 
there exists at least one constrained model o J ,i J and one pair of boundary conditions 
[Ti,72]( x>y ) such that Cj\^ 2 J x ^.V > 1. This means that there exists a constrained model, 
the one corresponding to a' , r', which does not fulfill the Dobrushin-Shlosman uniqueness 
condition DSU(V, 5) with 5 < 1. 

• / = 3: we have considered all possible constrained models and for each model we have 
considered 100 different pairs of boundary conditions. The results in Table 1 refer to a 
run of 1.3 • 10 6 sweeps. In this case there exists a particular constrained model a', r' and 
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a pair of boundary conditions [Ti,7^]( x y ) such that Up y, ^./ ( ^'f ) > 1 and Cpy^J^ < 1; 
that is, the upper bound is "too large", while the lower bound is "too low". Hence, for 
this model we can neither say that the Dobrushin-Shlosman condition is fulfilled nor the 
opposite; we must necessarily consider larger volumes. 

• I = 4: we have considered 50 constrained model and 60 pairs of boundary conditions; the 
results in Table 1 have been obtained by performing 1.3 • 10 6 full sweeps of the algorithm 
described in Section 4. 

• / = 8: we have performed the same statistics as in the case I = 4, but in this case it is 
obviously less significant, because the global number of possible choices is much greater. 
We have performed 1.3 • 10 6 full sweeps of our Monte Carlo. 

• / = 16: we have considered 30 constrained model and 30 pairs of boundary conditions; 
we had to reduce the number of runs, because of their length. 



Remark. 

The error AUpy in Table 1 is the usual statistical error, that is the standard deviation 
on the measure of the average Upy of the Monte Carlo results. On the other hand the 
best estimate Cpy and the error ACpy have been obtained by fitting the Monte Carlo 
results with the equation 

Cf>,v + —jj - , 

where I is the number of full sweeps of the run (see [BMO]) and ACpy = — -4^- . 



The results in Table 1 suggest that in the case / = 4 all constrained models satisfy 
the Dobrushin-Shlosman condition, that is the volume V with / = 4 is "large enough" 
for our purposes. But, strictly speaking, we cannot be sure about that, because we had 
to perform a statistics on the constrained models and on the boundary conditions; that 
is, there could exist a particular "bad" constrained model not satisfying the Dobrushin- 
Shlosman condition. Then we have considered values of / larger than 4 and we have shown 
(see Table 1) that for / = 8, 16 the value of Upy is so small and the effect of the change 
of conditioning spin is so localized (as we will explain) to lead one to the conclusion that 
the existence of such a bad model can reasonably be excluded. 

Our first observation refers to how the quantity jUpy, namely the average distance 
of the two copies of the system which evolve following the joint Monte Carlo dynamics, 
behaves as a function of /. Indeed, it grows from 0.530 to 0.8415 when one increases / from 
1 to 3, and then it remains approximately constant when one increases the value of /; this 
is what one expects heuristically. 
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In all the above described cases the "statistics" has been performed by choosing 
completely at random the constrained models and the boundary conditions; in all cases we 
have found that the most "dangerous" models, that is the constrained models with highest 
values of our upper estimator, are the striped and the chessboard ones. 

In order to strengthen the claim that inequality (2.23) is satisfied provided one 
chooses / large enough, in the case / = 6 we have performed a "rational" statistics, that is 
we have chosen the constrained models and the boundary conditions following reasonable 
criteria. We have chosen / = 6 for our final calculation, because the results listed in Table 
1 suggest that in this case condition (2.23) should be fulfilled, while, on the other hand, 
a full sweep of the Monte Carlo algorithm takes an acceptable CPU time so that we can 
perform a reasonably wide statistics. 

The algorithm introduced in Section 4 describes the "coupled evolution" of two copies 
of the same model, characterized by two different boundary conditions; we recall that at 
each instant of time t one and only one site (x, y) G V is updated and we observe that 
the property (4.5) suggests that the differences between the two copies of the model have 
a unique origin: the difference of the two boundary conditions in (x,y) G dV' + . 

Due to this fact it seems resonable to assume that during the evolution the total num- 



of the system is almost always equal to zero and sometimes these disagreements propagate 
in V starting from the site (x,y) in dV' + . Now, if one recalls that the average distance 
between the two copies is approximately equal to 0.87 for I > 4, it looks likely that the dis- 
agreements between the two copies of the system are localized in a "small" region around 



We have tested this hypothesis as indicated in the histogram in Fig. 4 which shows 
the spatial dependence of disagreements between the two copies. We have plotted the 
histogram for various constrained models, moving the site (x, y) along one of the four sides 
of V and for many pairs of boundary conditions; in all cases that we have considered, we 
have obtained histograms similar to that depicted in Fig. 4. The results summarized in 
Fig. 4 strongly suggest that the disagreements between the two copies of the model are 
almost completely localized in a suitably chosen 3x2 rectangular block R( x ,y) C V . 

Now, given the constrained model corresponding to a' G Q! v , and r' G Q' dv , + , the 
above remarks suggest that the average number of disagreements between the two copies of 
the system strongly depends on the values of a'^ x ^ with (x, y) G R( x ,y) an d weakly on the 
values of a'^ x y ^ outside R( x ,y)- Hence, we have performed the statistics on the constrained 
model in the case / = 6 by considering all possible constrained models only inside R( x ,y)- 
In the following we precisely describe how the statistics has been performed. 
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(x,y). 



• We have considered (x,y) = (0,3), R( x ,y) = {( x iV) £ V : x = 1,2 and 2 < y < 4}, 
we have modified the boundary condition in (x,y) in the six ways depicted in (5.1) and in 
each case we have considered two possible boundary conditions in dV' + \ {(x, y)}. 

• All possible constrained models have been considered in R^^, while in V'\R^^ we have 
considered only the chessboard model and the model with a', s = +1 W(x, y) G V'\R^,y)- 
Indeed, we expect that these two models are respectively the most and the less "dangerous" 
ones, as the results of the previous statistics suggest. 

• In each run of the joint Monte Carlo algorithm we have performed 10 5 full sweeps, that 
is we have updated the whole lattice V 10 5 times. 

The results can be summarized by saying that the most dangerous constrained models 
inside R( x ,y) appear to be the following ones 

+ -+ - + - +-+ - + -. ( , 2) 

+ - + - + - - + - +-+' (5 - 2j 

in particular, our (indeed quite small) statistics on the boundary conditions suggests that 
the most dangerous model among those in (5.2) is the second one; in this case, taking 
the chessboard model in V \ R( x ,y) we obtain Up y, 2 J X T V = 0.610. The numerical results 
confirm the weak dependence of the estimators on the constrained model outside R(x,y)i 
as well; actually the differences are of 5-10%. Finally, this set of Monte Carlo runs shows 
that the most dangerous ways in which one can modify the boundary conditions in (x, y) 
are the following ones 

+ +-- and +--+. (5.3) 

Once we have understood what are the worst situations inside R(x,y) an d in (x, y) we 
have performed the wide statistics on the possible boundary conditions described below. 

• We have considered the most dangerous constrained model inside R(x,y)- 

• We have considered ten possible constrained model outside R( x ,y)'i six of them are those 
depicted in Fig. 1, the remaining four have been chosen at random. 

• In each case we have considered twenty possible pairs of boundary conditions with the 
original spins in -S(x,y) and adjacent to V chosen like in (5.3). 

The weak dependence of the estimator on the constrained models outside R( x ,y) has been 
confirmed and we have found 

Upy, = 0.633 MApy, = 0.011 ; (5.4) 

hence we can confidently say that the condition (2.23) is satisfied if one considers the 
volume V with I = 6. 



1/ febbraio/2008 [25] 



5:5 



6. Conclusions. 

As we explained in Section 2 the problem of proving Gibbsianness of our renormalized 
measure is reduced to the verification of DSU(V, S) condition for some V and 5 < 1 for all 
possible constrained models. 

It is clear that disproving the condition for a given volume V is much easier than 
proving it since, to disprove, it is sufficient to exhibit one constrained model and one 
boundary condition for which a lower bound C for the uniqueness parameter appearing 
in DSU exceeds one; moreover since this lower bound involves variation distance and 
then thermal averages it is naturally computable via a Monte Carlo procedure. On the 
other hand an upper bound has to involve the consideration of all possible constrained 
models as well as all possible boundary conditions. Moreover a priori it was not clear 
how to provide an upper bound based on a Monte Carlo computation; this motivated 
the idea of the dynamical surgery. The necessity of a Monte Carlo approach comes from 
the consideration of how fast the number of constrained models and possible boundary 
conditions grows as a function of /, the side of the squared volume where we try to verify 
DSU. 

Since for / = 2 we find a particular constrained model and boundary condition for 
which C is greater than one, necessarily we have to go at least to / = 3 and already the 
number of independent computations is very large. Moreover since for / = 3 the lower 
bound seems always less than one whereas for at least one case the upper bound is larger 
than one we can neither disprove nor try to prove with our bounds the validity of DSU 
for a square with side 3 so that we have to go at least to I = 4. 

It appears clear from our numerical computations that our upper estimator U for 
the uniqueness parameter has the correct behaviour with /: as a consequence of the spatial 
localization of the set of disagreements between our two coupled processes U is inversely 
proportional to /; thus increasing I is the correct choice to get a value of U so smaller 
than one (included the error) to be confident in the validity of DSU . Unfortunately 
increasing / implies an enormous increase in the number of computations; introducing 
some statistics becomes necessary. The right compromise between smallness of U and 
number of constrained models and boundary conditions came out to be I = 6. In this case 
we performed the "rational statistics" that we have described in Section 5 by exploiting 
the (numerically evident) small dependence of U on the constrained model and on the 
boundary condition far from conditioning spin that we are changing. 

Our computations relative to the case / = 6 make us really confident on the uniform 
validity of DSU. 

We can say that our method is successful because even the correlation length of 
the "worst" constrained models is very small; however it is not so small to avoid the 
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consideration of sides / at least greater than 4. 

We want to make now some general remarks on our Monte Carlo algorithm. 

As it has been underlined in [DS1], the nice feature of finite size conditions like DSU, 
involving the behavior of Gibbs measures in finite volumes, is that they are able to imply 
absence of phase transitions and many nice properties of the unique infinite volume Gibbs 
state. This point of view is very helpful for example when we have to decide whether or 
not a given system is in the unique phase regime; especially when we do not have a natural 
parameter (like the inverse temperature (3) whose smallness imply weak-coupling. 

The "uniqueness test", based on the verification of DSU, has the advantage, w.r.t 
the traditional Monte Carlo test, of being based on rigorous grounds. However a real 
computer assisted proof seems very difficult to achieve unless the concerned volumes are 
really very small. If this is not the case the use of a Monte Carlo algorithm becomes 
essential; then the situation is somehow intermediate between a traditional Monte Carlo 
simulation and a computer assisted proof. 

We want to remark that our algorithm to compute a numerical upper bound on the 
Vasserstein distance between Gibbs measures, which is based on "local readjustment" of 
the coupling, seems to be quite performant and probably it can be used in more general 
contexts. Finally it is remarkable that, due to the very nature of the coupling procedure, 
the statistical error on U is much smaller that the corresponding one on C; indeed in this 
last case all sites of our volume and not only the disagreements, like in the computation 
of U, play a role as a source of statistical errors. 
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7. Appendix A. 

We consider the space S := {l,...,n} with n > 2 and the metrics p(s,s') = p(s,s') := 
1 — 5 SjS > Vs, s' G <§; let us denote by A and \x two probability measures on S and by K 
the set of the joint representations of A and \i. Hence, given q e /C one has that is a 
probability measure on S x 5" such that 



We recall, now, that the total variation distance and the Vasserstein distance between 
A and \i are respectively given by 





ses 



s'es 



Var(X, /i) 



-J>( S ) -/*(«)! 



(A.2) 




Proposition Al. 



1. With the notation introduced before 



V p (X,fi) = Var(X,n) 



(A3) 



2. Let us consider the following partition of the set S 



S=AUBUC , 



(AA) 



where 



A:={seS: A(s) >^(s)} 
B:={sGS: A(s) < fj,(s)} 
C:={seS: X(s) = fi(s)} 



(A.5) 



one has that 





seA 



Proof. 

1. See [Dl], page 472. 

2. It is an immediate consequence of the normalization of fx and A. □ 
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We want, now, to calculate the particular joint representation q* G K. such that the 
following equality is satisfied 



v p (x^) = p( s ' s ') ■ #*(«>«') ; 

s ,s'eS 



(A.7) 



that is we are looking for the joint representation "optimizing" the sum in the definition 
of V p (X,fj,). In other words, by virtue of the above Proposition, we can say that our aim 
is to find an n x n square matrix 



such that 



/ #l,n #2,n #3,r 



#1,3 #2,3 #3,3 
#1,2 #2,2 #3,2 
#1,1 #2,1 #3,1 



#n, n \ 



#n,3 
#n,2 

in, J 



(A.8) 



^2q* s ,s' = Vs> Vs' = 1, ...,n and ^ 9 * s , = A s Vs = 1, 



s'=l 



and such that the sum of the off-diagonal elements is given by 

s,s': s/s' s£A s£B 



?1 



(A9) 



(A10) 



we have introduced the notation: A s := A(s) and /i s := fx(s) Vs G 5". 

In order to compute all the elements of the matrix g* s , one can follow the procedure 
described below; first of all we introduce the following filling rule. It is a construction of 
a matrix q* along a sequence of steps. At each step some new set of entries of q* will be 
determined whereas for the rest of the entries only some suitable sums of them will be 
determined. We will end up with a matrix q* satisfying both (A. 9), (A. 10) and in this way 
we show that our procedure is well defined in the sense that at each step it is not empty 
the set of matrices satisfying our requirements. 

1. If A s = n s Vs G S then 

#s, s ' = \s8.s,s> Vs, s' G S 

2. If 3s G S such that A s ^ fj, s then Vs G S one puts 

# s *, s = min{A s ,^ s } 
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and 



if fi s = X s then 

if n s > X s then q* , 




if A s > fi s then q*, 8 



s': s'^s 

Vs' 7^ s and g*, s are subject to the condition : 



s': s'^s 



In Fig. 5 an example of application of the above filling rule is shown. 

In case 1 the problem is trivially solved; indeed by means of the above filling rule 
one obtains a diagonal matrix fulfilling the requirements (A. 9) and (A. 10). In case 2 it 
is immediate to see that we construct a matrix satisfying the requirements (A.9), but the 
problem is that in general not all off-diagonal terms are determined: many of them are set 
equal to zero, but for the others only some sums are specified and not the single values. 

By collecting all off-diagonal no n- vanishing entries one obtains an |^4| x \ B\ rectangu- 
lar matrix, whose elements are not yet determined. More correctly we should speak about 
a class of matrices (which at the end of our construction will be shown to be non-empty): 
the ones satisfying the requirements corresponding to one step of application of our filling 
rule. 

Let us denote by q^j with 1 < i < ri\ := \A\ and 1 < j < n[ := \B\ the elements of 

this new matrix. We observe that qf- corresponds to the element q^^ of the original 
matrix, where 



It is important to observe that, due to the filling rule, the elements of the new matrix 
must satisfy the following requirements 



/(I): 
/(*) = 



min{s E S : X s > fx s } 

min{s > f(i — 1) : A s > /i s } Vz > 1 



(All) 



and 



f'd) ■ 



min{s E S : jj> s > X s } 

min{s > f'(i — 1) : jJL s > X s } \/i > 1 



(A12) 
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From the fact that has been constructed by collecting all off-diagonal elements of q* 
different from zero and from equations (A. 13), it follows that the (A. 10) is fulfilled. 

In order to fill the matrix q^ (by further reducing its indetermination) one applies 
again the filling rule stated above. Again two things may happen: all elements of q^ are 
fixed or one obtains another 77-2 x n' 2 rectangular matrix q^ to be filled. 

One goes on applying the filling rule and in this way one obtains a sequence of 
rectangular matrices 

q*, q W, q W,..., <?«... 
which are smaller and smaller in the sense that 

n ■ n > ri\ • n! x > 712 ■ n' 2 > ... ; 

the procedure is stopped when all elements are fixed. 
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Table 1 



/ 








ACpy 


2 


1.452 


0.0074 


1.34 


-2.50 


3 


1.122 


0.0071 


0.90 


+2.05 


4 


0.877 


0.0046 


0.73 


+0.32 


8 


0.436 


0.0026 


0.30 


+0.74 


16 


0.207 


0.0016 


0.16 


+0.04 



Table 1: In the table we have listed the values of Upy , AUpy, 
Cpy and AjCpy obtained at I = 2,3,4,8,16. The meaning of 
AUpyi and ACpy is explained in the Remark in Section 5. 
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Fig.l: Examples of configurations of renormalized variables corresponding 
to constrained models which are particularly important in our calculations. 
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Fig.2: Specific heat for the chessboard model defined in (3.1) as a function 
of the inverse temperature (3. 
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Fig.3: Specific heat for the striped model defined in (3.1) as a function of 
the inverse temperature (3. 
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Fig. 4: Each square of the plane X— Y represents a site (x, y) of a lattice V with / = 6; the X 
axis is oriented from the left to the right, while the Y axis is oriented from the top of the figure 
to the bottom. The results contained in the histogram refer to the chessboard model: the height 
of each column (x, y) is given by the ratio 

(n,v) (x,y) 

El=i E(x, y)e i/'( 1 - ^w.;^) 

where 1=1- 10 5 is the number of full sweeps of the run, S^}'^ and are the block-variables 

of the two copies of the system on the site (x, y) and after s sweeps. The results in the histogram 
have been obtained considering a pair of boundary conditions [%_, 7^] (0,3) j such that the two pairs 
of original spins in the block -B(o,3) are H — and — h 
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A 4 
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M3 
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A 2 
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n 
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Ai 


A 2 


A 3 


A 4 


A 5 



a (1) n W a {1) n* \a 

<?1,2 12,2 %,2 M4 - M 

Ql,l 12,1 %,1 M2 - A 2 



Ai — [i\ A3 — fis Ag 



M5 



(a) 



(6) 



Fig. 5: Application of the filling rule in the case: n = 5, Ai > fi\, A 2 < /U 2 , A3 > /X3, A4 < ^4 
and A5 > /U5. In (a) we have depicted the matrix q* s ,: the sum of the elements of one of its rows 
must be equal to the corresponding value depicted on the left, while the sum of the elements of 
one of its column must be equal to the corresponding value depicted on the bottom. The squares 
represent elements which are not fixed after the filling rule is applied for the first time. The 
collection of all squares gives rise to the matrix q^j which is depicted in (b). The values in the 
right-hand column and in the bottom row are the constraints (A. 13) which must be satisfied by 
the elements, respectively, of each row and each column of the matrix. 
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